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INTRODUCTION 


The following tables contain the neutral temperature, neutral 
densities for N 2 , 0 2 , 0, Ar, He and H, mean molecular weight, and total 
mass density as predicted by the MSIS (Mass Spectrometer and Incoherent 
Scatter) empirical thermosphere model [Hedin et al . , 1977a; Hedin et al., 
1977b] for selected altitudes, latitudes, local times, days and other 
geophysical conditions. The model is based on a least squares fit to 
density data from mass spectrometers on five satellites and temperature 
data from four incoherent scatter stations, providing coverage for most 
of solar sunspot cycle 20. 

Included in the model data base are longitudinally averaged N^, He, 
and 0 densities from the OGO-6 mass spectrometer [Hedin et al., 1974], 
longitudinally averaged N 2 , He, 0 and Ar densities from the AEROS-A 
(NATE) mass spectrometer [Spencer et al., 1974], the N 2 , He, 0, and Ar 
densities from the San Marco 3 mass spectrometer [Newton et al., 1974; 
Newton et al., 1975], the N 2 densities from the AE-B mass spectrometer 
[Reber et al., 1971], and the N 2 , He, 0, and Ar densities from the 
AE-C (OSS, NACE, NATE) mass spectrometers [Nier et al., 1973; Pelz et 
al., 1973; Spencer et al., 1973]. The 0 2 and H densities are inferred 
using ion mass spectrometer data from AE-C (BIMS) [Brinton et al., 

1973] . Neutral exospheric temperature data are included from Arecibo 
[Waldteufel and Cogger, 1971], St. Santin [Alcayde, 1974], Millstone 
Hill [Salah and Evans, 1973], and Jicamarca [McClure, 1971]. 

The total number of data points available (over 100,000 for some 
constituents) is much larger than needed to define the model parameters 
Thus a subset of approximately 2500 points per constituent was selected 



by sorting the data into "boxes" based on altitude, latitude, local time, 
day of year, mean 7 , 27 day 7 variations, and magnetic activity, 
and then selecting data at random irom these boxes. The purpose of this 
procedure is to spread the selected data points as evenly as possible 
over the relevant geophysical parameters. No averaging is performed. 

The model formulation is based on the Jacchia [19651 formula for 
temperature profiles and the related diffusive equilibrium density 
formula [Walker, 1965]. The exospheric temperature T — and other quantities 
are represented by expansions (G(L)) in terms of spherical harmonics in 
local time and latitude plus additional terms for variations with the 
day of year, solar flux index Fjq 7 , and magnetic activity index Ap. 

T(Z,L) = TJL) - [TJL) - T 120 (L)1 exp[-s(L) (Z-120)l 
n(Z,L) = n 120 (L) DfZ.TJL), s(L)l 
D[Z, TJL), s(L)l = [(l-a)/(l-ae"° 5 )l 1+ ° +Y exp(-o Y C) 

where. 


TJL) = TG(L) 

o = s(L) + .00015 

T 120 (L) S f 120 G(L) 

£ = (Z-^OHRg+^Ol/CRg+Z) 

s(L) - ?G(L) 

Y = Mg^/fRTJDol 

n 120 (L) " "l20 •*pl®C L 3-ll 

a - 1 -T 120 (L) /T^CL) 

G(L) ■ formula for the variation of temperature and other quantities 

as a function of local time. 

latitude, day of year, F 1Q 7 , 


Ap and numerical coefficients determined in the modeling 
process. 


iv 


g J20 = 9.44663 x 10' 3 km/s 2 

R_ = 6356.77 km 
E 

R = 8.314 x 10" 3 g km 2 /mol s 2 K 

L = argument list of geophysical and geographical parameters (local 
time, latitude, day of year, and Ap) 

Z = altitude, km 

a = thermal diffusion coefficient, equal to -.4 for He and H and 
0.0 otherwise 

M = mass number of gas species 
T = ambient temperature, K 
n = ambient density, cm" 3 
s « temperature gradient parameter, km~* 

Ti 2 q = average 120 km temperature, K 
T = average exospheric temperature, K 
s = average gradient parameter, km"* 

”l20 = avera 8 e density at 120 km, cm" 3 

The diurnal and magnetic activity variations are multiplied by 
additional height dependent terms to allow for departures from diffusive 
equilibrium below 200 km. 

The 120 km temperature coefficients are specified first on the basis 
of a very limited amount of published and unpublished incoherent scatter 
data. Next the 120 km N 2 density coefficients are specified by fitting 
low altitude data from AE-C and San Marco. The exospheric temperature 
and the temperature gradient coefficients are determined next by fitting 
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the exospheric temperature and ^ density data taken on days with Ap 
less than 10 and then determining the magnetic activity related parameters 
from data at all levels of magnetic activity (Ap = 2 to 130). The model 
coefficients for He, 0, and Ar which specify the effective density 
variations at 120 km are determined using data from above 190 km, but 
departures from diffusive equilibrium are based on data down to 140 km. 

The 0 7 and H coefficients are determined from data in the 250 to 400 km 
altitude range. 

The MSIS model is most representative of magnetically quiet conditions 
with best data coverage between 200 and 600 km and with solar Fj 0 7 
indices between 75 and 180. While only one third of the data are associ- 
ated with an Ap greater than 10 and no attempt was made to predict short 
term changes, the predictions for high Ap should be representative of 
the conditions to be expected during magnetic storms. Model densities 
should be reasonably valid above 150 km, but temperatures below 200 km 
may not adequately represent the diurnal and semidiurnal variations. 

The model densities values below 150 km are essentially extrapolations 
from higher altitudes. The data coverage and quality of fit for the 
various gas species is far from uniform. A crude ranking of overall 
confidence based on coverage, quality of fit, and number of data points 
is as follows: Temperature and Nj* 0, He, Ar, H, and C^. The absolute 

error in the mean temperature and mean density profiles is approximately 
+2, *15, +20, +20, +25, +50, +50% respectively. Hydrogen densities 
below 200 km are expected to depart strongly from diffusive equilibrium 
and this is not reflected in the tables. Both H and Oj variations are 
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based on inferred densities (from ion composition) and need confirmation 
by direct measurements as they become available. 


The Tables have been generated for five latitudes (90°, 45°, 

0°, -45°, -90°), four local times (3, 9, IS, 21 hours), four days of 
the year (81: equinox, 173: summer, 265: fall, 356: winter), four solar 
F 10 7 levels (75, 100, 150, 200), and three levels of Ap (4, 20, 100). 
The table headings are defined below: 

LAT = Geographic latitude (deg) 

LT = Local solar times (hrs) 

DAY = Number of days since January 0 

F107 * Mean 10.7 radio flux (10" 22 Wm" 2 H "*) over three 

Z 

solar rotations centered on period in question*. 

Ap * Daily magnetic index for day in question 

ALT = Altitude (km) 

TEMP * Neutral temperature (K) 

MW • Mean Molecular Weight 

MASS DENS = Total Mass Density (gm/cm 3 ) 


*The tables do not explicitly allow for the 27-day solar rotation 
effect but this can be approximated by using F107«.75 Fjq 7 + .25 Fjg 7 
where Fjq 7 is the mean flux and F^ ? is the daily flux for the 
previous day. 
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